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The silicon vacancy is a prominent intrinsic defect of cubic SiC (3C-SiC) to which much effort 
has been devoted so far, experimentally and theoretically. We calculate its properties using the 
GW approximation that does not suffer from the band gap problem. The obtained formation and 
transition energies deviate significantly from the usual density functional theory evaluations and 
now compare favorably with experiment. A new assignment for the main line of photoluminescence 
is then proposed. We further perform GW calculations for the saddle point of reaction paths. The 
resulting barrier energies explain the thermal annealing experiments thanks to an original mechanism 
mediated by a minority charge configuration. 



I. INTRODUCTION 

Cubic silicon carbide (3C-SiC) is considered as a 
promising semiconductor for use in severe environ- 
ments. The potential applications include hardened 
electronic devices and nuclear fuel coatings. In par- 
ticular, its temperature and radiation resistance makes 
3C-SiC attractive for the next generation of high- 
temperature gas-cooled nuclear reactor^. Due to these 
purposes, irradiation-induced point defects of 3C-SiC 
have been widely studied both experimentally^ and 
theoreticall}'^^. Among the intrinsic defects, the silicon 
vacancy plays a prominent role since it remains stable 
up to relatively high temperatures and it is easily identi- 
fied by electron paramagnetic resonance (EPR) measure- 
ments due to its high-spin configuration^ . Unfortunately, 
a careful comparison between the experimental data and 
the ab initio calculations obtained within density func- 
tional theory (DFT) is not satisfactory. 

The usual approximations of DFT suffer from the in- 
famous band gap problem. There is nowadays a com- 
mon agreement that the band gap problem particularly 
plagues DFT predictions of defect properties in semicon- 
ductors and insulators^. Point defects in non metallic 
solids generally produce additional electronic levels in- 
side the band gap region. The precise location of these 
levels is critical for most defect properties: charge state 
transition energy, formation energy, photoluminescence 
(PL) lines. Predictive calculations absolutely need to be 
based on schemes that are devoid of the band gap prob- 
lem. The GW approximation of the many-body pertur- 
bation theorjEI has be en pr oven for several years to yield 
the correct band gap a^°^ . Still the GW calculation of 
defects is a great numerical challenge. 

In the present article, we calculate, using the GW ap- 
proximation, the energetics and the stability of the silicon 
vacancy in 3C-SiC. In order to access these properties, 
we performed the GW calculations not only for equilib- 
rium structures, but also for saddle point positions. The 
GW results were combined so to reduce an error, named 
concavity, that was recently identified for the GW ap- 
proximation by one of uJ^. The procedure, explained in 
Sec. |lll allows us to calculate the formation energies of 



the metastable form Vg; and of the stable configuration 
(VcCsi) within the GW framework in Sec. Ill The re- 



sults bring a quantitative agreement with various exper- 
imental measurements: EPR, PL, and thermal anneal- 
ing. Noticeably, the calculations propose a reassignment 



of the main PL line in Sec. IV and highlight an origi- 
nal decay mechanism through a minority charge state in 
Sec.lVl 



II. METHOD 
A. GW approximation applied to the defects 

Thanks to its accuracy in predicting band gaps and 
electronic levels, the GW approximation appears to be 
the method of choice to deal with point defects in semi- 
conductors and insulators. However some noticeable bot- 
tlenecks have inhibited the use of this scheme so far. Very 
recent workJi^ES! have lifted the theoretical and numeri- 
cal problems as recapitulated in the following. 

Deriving from Green's function theory, the GW ap- 
proximation provides meaningful quasiparticle electronic 
levels. The GW approximation is hence reliable in pre- 
dicting charge changes with constant geometry. For in- 
stance, the energy of the highest occupied molecular or- 
bital (HOMO) can be interpreted as a total energy dif- 
ference at constant geometry, which in turn defines the 
ionization potential /: 



(1) 
(2) 



where the first argument refers to the geometry of the 
defect and the second shows the actual charge state q. 
Eq stands for the ground-state total energy and e for 
the quasiparticle energy. Symmetrically, the lowest un- 
occupied molecular orbital (LUMO) energy is a total en- 
ergy difference, which can be understood as the opposite 
of the electron afhnity: —A. In the context of defects, 
these quasiparticle energies are named vertical transition 
energies ey{q/q + 1). Note that this vertical transition 
energy ev{q/q + 1) could be obtained in principle either 
from eg^MO ( ^si > ?) or from e^^^ ( V|i ,q + l). One of us 
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showed recentljff^l that these two values slightly differ in 
practice due to the concavity of the GW approximation. 
The most consistent way to handle this discrepancy is 
hence to approximate the vertical transition energy by 
the mean value of the two mentioned quasiparticle ener- 
gies. 

The calculation of defect formation energies also re- 
quires to relax the geometries, since the structure of a 
defect changes according to the charge state. Unfortu- 
nately the GW total energies and forces are out of reach 
even in state-of-the-art implementations. That is why 
Rinke and coworkersP^l recently introduced a scheme that 
combines DFT and GW. DFT is used to deal with the 
structural changes at constant charge state. The GW ap- 
proximation is employed for charge changes at constant 
geometry. When traveling on a given Born-Oppenheimer 
surface, the band gap underestimation is expected to be 
harmless and therefore DFT can be safely used. Then 
the GW approximation allows one to address the ver- 
tical transitions from one Born-Oppenheimer surface to 
another one. The combination of the two moves finally 
permits one to change both the geometry and the charge 
states. The thermodynamic charge transition energies 
^thil/Q + 1) can thus be obtained by introducing the en- 
ergy of an intermediate point: 

eth{q/q + 1) = Eo{Yl„ q) - Eo{Yl+\q) 

+Eo{Vlt\q)-EoiYlt\q+l)- (3) 

The first two terms account for a structural change at 
constant charge q, which is obtained accurately from two 
DFT calculations. The last two terms account for a 
charge change at constant structure, which is calculated 
from GW quasiparticle energies as explained above. 

With the present combined DFT /GW approach, the 
absolute formation energy of the defect for any charge 
state can be obtained, provided one reference formation 
energy. In Ref. |T2| the reference formation energy was 
chosen to be the DFT total energy with all defect states 
in the band gap empty, so that the band gap problem 
does not enter in the total energy. Then the neighboring 
charge state formation energy can be obtained as exem- 
plified here: 

EfiYti) = EfiVlt) + etkil + /2+) - e^g^i - Me, (4) 

where Ef stands for the usual formation energy and fi,, 
stands for the Fermi level with its zero set at the GW 
valence band maximum, eyBM- 



B. Technical details 

The main numerical bottleneck of the GW approxima- 
tion comes from the dependence of the GW exchange- 
correlation self-energy onto the empty states. This prob- 
lem was given a partial answer in Ref. 1151 which allows 
us to run GW calculations for supercells as large as 215 
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FIG. 1. (Color online) Formation energy (silicon-rich) of the 
siHcon vacancy in the charge state 2+ as a function of the 
supercell size. We provide the uncorrected data with open 
square symbols, the potential aligned data with triangle sym- 
bols, the Madelung monopole corrected data with circles, and 
the Lany-Zunger corrected data with diamonds.^* In addition, 
the dashed lines represent reasonable extrapolations based on 
the three largest supercells. 



atoms using a very low 3:1 ratio between empty and oc- 
cupied states. 

We employ the usual perturbative GW method 
{GoWoY^* with norm-conserving pseudopotentials and 
the Godby-Needs plasmon-pole model, using a 30 Ha cut- 
off for the wavefunctions and a 6 Ha cutoff for the dielec- 
tric matrix. The k-point grid is a 2x2x2 Monkhorst-Pack 
grid for the DFT force calculations. It is reduced to the 
F point only when turning to GW calculations. The cal- 
culations are spin polarized and the neutral vacancy is 
approximated by its triplet spin state, which m ay oye res- 
timate its formation energy by less than 0.1 eV.C^The 
obtained band gap is 1.35 eV within the local density ap- 
proximation (LDA), 2.19 eV within GW to be compared 
to the experimental value, 2.37 eV. The discrepancy be- 
tween the GW band gap and the experimental band gap 
can be considered as the error bar of the results presented 
in the following. 

The calculations of charged defects in replicated su- 
percells notoriously converge slowly as a function of the 
supercell size. The supercell technique gives rise to im- 
portant finite-size effects, among which the electrostatic 
interaction between the charged defects is arguably the 
largest. In order to address this issue, a Madelung cor- 
rection can be designed.— However it has been shown 
in the recent years that a straightfor ward M adelung cor- 
rection does more harm than benefit.E^I^SI^ We clarified 
this issue for our specific system by performing a careful 
convergence study about the silicon vacancy Vg^, which 
should be the most dramatic case of our study. 

In Fig. [ij the formation energy of Ygf is provided for 
different cubic supercells ranging from 63 atoms to 999 
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FIG. 2. (Color online) Formation energy of the silicon vacancy 
Vsi (top) and of the complex (VcCsi) (bottom) in 3C-SiC 
under silicon rich conditions within LDA (dashed lines) and 
within the GW approximation (solid lines) as a function of 
the Fermi level /Xe- The silicon atoms are light gray (yellow) 
and the carbon atoms are dark gray (blue). 



atoms. The structure is frozen with all the atoms in 
the perfect crystal positions except for the four nearest 
neighbors to the vacancy, whose positions were relaxed 
in a 63 atoms supercell. This explains why the formation 
energy is slightly larger than the fully relaxed energy re- 
ported elsewhere in the text. In the present case, the 
Madelung correction clearly overestimate the correction. 
In Fig. [l] we provide the curve corrected by the popu- 
lar Lany-Zunger correction schem^^^, which consists in a 
potential alignement together with a reduced Madelung 
correction (precisely 2/3 of the usual Madelung term). 
The Lany-Zunger scheme is also not satisfactory as shown 
by the figure. Note that the dielectric constant has been 
calculated ab initio consistently using the same param- 
eters as the defect calculations. We then found more 
reliable to perform a mere potential alignment. Using a 
215 atoms supercell with no charge correction allows us to 
evaluate within 0.1 eV the formation energy ofV|+ The 
values compared to experiments are always differences of 
formation energies, which are expected to be even more 
accurate, since the errors compensate to some extent. 



III. FORMATION ENERGY OF THE 
METASTABLE SILICON VACANCY AND THE 
STABLE COMPLEX 

The silicon vacancy is a high energy defect that can 
only be observed in heavily irradiated SiC. The numerous 
silicon vacancies created by irradiation are a metastable 
configuration: the silicon vacancy experiences a large en- 
ergy drop when transforming into the complex made of 
a carbon vacancy and a carbon antisite (VcCsi)^. Fur- 
thermore, the experimentally observed silicon vacancies 



possess a pecuHar high-spin configuration (S=3/2) with 
three aligned spins resulting in a 1- charge state. This 
particular spin configuration gives rise to an unambigu- 
ous EPR signal, named Tl center. The Tl center can 
be observed in irradiated samples, whatever the initial 
doping conditions and irradiating particle^ . After ir- 
radiation, due to the quantity and the variety of point 
defects introduced in the material, the donor and accep- 
tor defects give rise to a compensated system and the 
Fermi level remains pinned in the vicinity of mid-gap. 

In Fig. [2] we provide the formation energy of the two 
competing configurations for the silicon vacancy: the 
original silicon vacancy Vg; and the complex (VcCsi). 
The LDA clearly suffers from the band gap problem so 
that all the transition energies are constrained to be be- 
low the LDA conduction edge, symbolized by a vertical 
dotted line. In addition, it would be doubtful to deter- 
mine the mid-gap region: would it be the middle of the 
Kohn-Sham band gap or the middle of the experimen- 
tal band gap? Conversely, the GW results offer a much 
clearer view on the system. First, the GW transition 
energies are correctly placed within the full range of the 
GW band gap, which matches reasonably well the exper- 
imental band gap. Second, the silicon vacancy adopts a 
correct 1- charge state in a high-spin configuration for a 
wide region around mid-gap. The silicon vacancy appears 
as metastable with respect to the complex (VcCsi) for 
any Fermi level in the band gap. The energy differences 
are sizable, ranging from 4.4 eV for the p-type region to 
0.7 eV for the n-type region. The GW results show that 
the complex carries a 2+ charge state in the major part 
of the band gap. 



IV. PROPOSAL FOR AN INTERPRETATION 
OF THE E-LINE OF PHOTOLUMINESCENCE 

The zero-phonon line of PL spectroscopy measures the 
vertical transition energy down from an excited defect, 
as schematically shown in the left panels of Fig. [3j In 
irradiated samples, additional features appear that can 
be linked to the intrinsic defects. In particular, a line, 
labeled E by Itoh and coworkerPI, was convincingly at- 
tributed to the silicon vacancy. Indeed the E line and the 
above mentioned Tl center in EPR have precisely the 
same behavior upon thermal annealing: the same three 
annealing stages. This clearly demonstrates that the pho- 
tons of the E line with an energy of 1.91 eV arise from 
electronic transitions occurring at the silicon vacancy. 

In Ref. 131 the authors proposed that the E line comes 
from the recombination of a conduction electron with a 
hole on a silicon vacancy, i.e., Vg;. So that they ten- 
tatively assigned it to the transition level et,(0/l— ) at 
an energy of i?c — 1-91 eV. This assignment was mainly 
supported by calculations: for instance, our LDA results 
place the transition eu(0/l— ) ^ Ey + 0.54 eV. By com- 
bining the LDA transition level together with the exper- 
imental band gap, 2.37 eV, the resulting emitted photon 
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FIG. 3. (Color online) Transitions associated to the photolu- 
minescence of the silicon vacancy in 3C-SiC. The schematic 
Born-Oppenheimer surfaces are represented on the left-hand 
panels. The PL levels calculated within the GW approxima- 
tion are represented on the right-hand panels. 



matches the experimentally measured photon. 

However, this assignment is questioned by our calcu- 
lations. When performing higher accuracy GW calcu- 
lations with the fair GW band gap value, the obtained 
vertical transition lies at e^(0/l— ) = + 0.79 eV = 
Ec — 1.58 eV, as represented in Fig. [3j In this improved 
framework, the obtained transition energy does not fit 
the experimental interpretation anymore. The disagree- 
ment is definitely larger than the calculation uncertain- 
ties, like the slight underestimation of the band gap by 
the GW approximation (2.19 eV instead of 2.37 eV) and 
the neglect of the excitonic effects. 

From our calculations, another interesting scenario can 
however be drawn that challenges the previous interpre- 
tation of the E line. In fact, the PL technique is not 
capable of distinguishing donor-valence recombinations 
from acceptor-conduction ones. So that the observed 
E line may rather arise from the recombination of an 
extra electron on the silicon vacancy (i.e., Vgj") with a 
hole in the valence bands. Then, the electronic transi- 
tion that emits the PL photon will be a signature for 
the vertical transition e^{2 — /I—). Previous evaluations 
for this transition were biased by the band gap prob- 
lem: within LDA, e„(2 - /!-) = Ey + 1.15 eV close to 
the LDA conduction edge (1.35 eV) and therefore was 
disregarded for the E line. Within GW, this transition 
is located interestingly close to the experimental E line. 
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FIG. 4. (Color online) Energy path of the transformation 
of Vsi into (VcCsi) for different charge states. The carbon- 
rich conditions were set. The Fermi level has been set to 
the experimental mid-gap value. The upper panel shows the 
5 image CI-NEB path obtained for the 215 atoms supercell 
within LDA. The lower panel provides the GW evaluation for 
the energy of the 3 critical points of the path. The points 
actually calculated are depicted with the symbols and the 
lines are guides to the eye. 



e„(2 — /I—) = Ey + 1.85 eV. A small error of the pho- 
ton energy is in fact expected, since our single electron 
picture does not take into account the exciton binding 
energy and since our calculated band gap is slightly too 
small. Note that the two sources of error have different 
signs. Indeed, the electron-hole attraction always low- 
ers the emitted photon energy. Exciton binding energies 
comprised between 0.18 and 0.24 eV were recently cal- 
culated for the carbon vacancy in 4H-Si CP. It is likely 
that an increase of the band gap would push the defect 
level up, since it is close the conduction edge. Our data 
clearly advocate for a reinterpretation of the nature of 
the E line: the E line is a signature of the transition 
from charge state 2- to charge state 1-, with an exper- 
imental value of £v{2 — /I—) = Ei, + 1.91 eV against a 
GW calculated value of e„(2 - /1-) = E^ + 1.85 eV. 



V. ANNEALING OF THE SILICON VACANCIES 
THROUGH A MINORITY CHARGE STATE 

Along with the low temperature EPR and PL studies, 
the thermal annealing o f the silicon vacancy signatures 
have been widely studiecplil22l The annealing of the sil- 
icon vacancy related signals occurs in three stages. The 
recovery stages I and II at rather low temperature (400- 
700 K) anneal about half of the signal, then stage III at 
high temperature (1050 K) marks the complete annihila- 
tion of Vgj. The activation energy of the last stage was 
determined to be Ea = 2.2 ± 0.3 e^^^. The thermal 
annealing curve yields accurate effective activation ener- 
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gies, but do not give any clue about the mechanisms at 
the atomic scale. In the following we propose a picture 
of the recovery stages of Vgj thanks to the published ab 
initio data and our GW results for a reaction path. 

It is reasonable to assume that the first two recovery 
stages are associated to the elimination of silicon vacan- 
cies by recombination with carbon and silicon intersti- 
tials, whose migrati on an d recombination barriers are in 
the range 0.5-1.5 eV'^HHl, "We have to exclude from this 
process the 4+ charged silicon interstitials, whose migra- 
tion energy is expected to be large. Vacancies themselves 
can also be considered immobile^. The extinction of the 
T5 center in EPR, assigned to Cgpu^'^, at the same tem- 
perature as the first stage suggests that this is related to 
the Vg; -I- Cgpjjj — >■ Csi recombination. The antisite is a 
very low energy defect that is very stable once formed. 

The third recovery stage shows a first order kinetics 
according to the isothermal annealing of Ref. H] This 
stage was tentatively proposed as a transformation of the 
metastable Vg; into the lower energy complex (VcCsi)'^^, 
but our calculated barriers do not agree with the exper- 
imental activation energy as shown in the following dis- 
cussion. Again we argue that the calculations suffer from 
the band gap problem and that the mechanism involved 
is more complex than was thought previously. 

We therefore performed reaction path calculations 
within LDA for the relevant charge states. The LDA 
configurations for the saddle points were then used to 
build up the GW formation energies as described above. 
The LDA paths were obtained thanks to the climbing 
image nuclg ed elastic band method (CI-NEB)2^ with 5 
images for 215 atoms supercells as displayed in the upper 
panel of Fig. [4j In order to compare the different charge 
states, we set the Fermi level to mid-gap. One can safely 
assume that the Fermi level is pinned by the numerous 
deep donors and acceptors, which somehow compensate 
each other in heavily irradiated materials. In focusing 
on the Fermi level we implicitly assume that the ther- 
modynamic equilibrium between charge states is always 
satisfied. This hypothesis is not true in general. However, 
in the present case, the annealing experiment is carried 
out at very high temperature (1050 K) and the kinetics 
of charge equilibration is a matter of nanoseconds. We 
will discuss this point in further details later on. The 
obtained LDA transition paths have a minimal barrier of 
2.9 eV for the 1- charge state. This value deviates much 
from the experimental activation energy. 

In the lower panel of Fig. |4] we provide the GW en- 
ergies for saddle points and for stable points. The GW 
barriers shine light onto a very peculiar diffusion mecha- 
nism. The vast majority of the silicon vacancies possesses 
a 1- charge state, however the direct transformation into 
the complex (VcCgi)" is blocked by a large barrier of 
2.75 eV. The very few positively charged silicon vacan- 
cies Vg; experience a low energy barrier of 1.50 eV. How- 
ever the total effective activation energy of the reaction 
in the 1+ charge state is the sum of the barrier energy for 
charge 1-f and of the energy needed to turn the vacancy 



from 1- to 1+. We finally obtain an activation energy of 
2.32 eV, which lies within the experimental uncertainty. 
Some lowering of this figure by temperature effects can be 
further expected. The transformation of the metastable 
silicon vacancy into the stable complex is mediated by 
a double charge change. This extinction of Vg; is inter- 
estingly driven by an ultra minor charge state: even at 
1050 K, the concentration of the 1-f charge state is only 
[Vsi]/[Vsi] 10"^- Note that only the GW values for 
the formation energies of Vg; and of the barrier heights 
were able to explain the experimental observations. 

Let us now assess the validity of instantaneous equi- 
librium between the different charge states of Vg;. All 
our calculated paths are based on this assumption. In- 
deed, this hypothesis can be validated by experimental 
data from Deep Level Transient Spectroscopy (DLTS) 
on a very similar system, hexagonal SiC (4H-SiC). 4H- 
SiC is expected to resemble much 3C-SiC, except for the 
band gap which is noticeably larger. DLTS experiments 
performed in heavily irradiated 4H-SiC measure a gigan- 
tic timescale to restore charge equilibrium of about two 
weeks for deep levels at room temperature.^^ These lev- 
els have been measured 1.2 eV away from the band edges 
of 4H-SiC, which constitutes an upper bound for a deep 
level in 3C-SiC, for which the band gap is only 2.37 eV. 

The theory underlying DLTS specifies that the recov- 
ery time T follows an exponential law*2Si 

T(xexp{AE/kBT), (5) 

where AE is the energy difference between the defect 
level and the closest band edge, fcs is the Boltzmann 
constant and T the temperature. In the annealing exper- 
iments of Itoh and coworker^, the critical temperature is 
1050 K. For this high temperature, the exponential law 
in Eq. ([5| gives a surprising collapse of the recovery time, 
which value is only 2 ns. The annealing experiment usu- 
ally lasts for several minutes or hours. As a consequence, 
the charge re-equilibration is many orders of magnitude 
faster than the annealing experiment itself. The instan- 
taneous charge equilibration can be safely assumed as we 
did to draw Fig. |4] 

VI. CONCLUSIONS 

As a conclusion, the use of the state-of-art GW ap- 
proximation for large supercells of 215 atoms allowed us 
to completely clarify the experimental observations con- 
cerning the silicon vacancy in 3C-SiC. The calculated 
properties deviate noticeably from the previous LDA 
studies, owing much to the absence of the band gap prob- 
lem. Beside the nice agreement with the experimental 
data, we brought deeper understanding in the physics of 
the silicon vacancy: we proposed a new assignment for 
the main PL line, which should be attributed to the 2-/1- 
charge transition; we confirmed that the third recovery 
stage of Vg; is indeed a transformation into (VcCsi) and 
proposed a mechanism that channels through the very 
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rare 1+ charge state. This last apphcation was possible 
only using GW calculations for the reaction saddle point. 
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